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1) Label root voxeJ with integer 0; 

2) Construct a queue and line up the root in the queue; 

3) If( There is at least one voxel in the queue ) 

Serve the voxel x on the top of the queue: 
For( each of x's 26-connected neighbor voxel y) { 
Iff y in the volume and has not been labeled yet ) { 
Line up the y in the queue; 

/* label the voxel y */ 
Set dist s= 999999 ; 

For( each of y's 26-connected neighbor voxel z ) { 

If( z in the volume and has been already labeled with an integer of n z ) { 
d z =n z +d(y,z); 

where d(y,z)is 10, 14 or 17 if the Euclidean distance between y and z is 
, respectively; 

} 

U(dist>d z ) { 

label y with integer dist ; 
dist = d z ; 

} 

J 

} 

) 

x leaves the queue; 

} 

Else { 

end of calculating the distance map 

} 



